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We consider stimulated Raman adiabatic passage (STIRAP) processes in tripod systems and show 
how to generate purely geometric phase changes of the quantum states involved. The geometric 
phases are controlled by three laser fields where pulse shapes, relative field strength and phases can 
be controlled. We present a robust set of universal gates for quantum computing based on these 
geometric phases: a one-qubit phase gate, a Hadamard gate and a two-qubit phase gate. 
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I. INTRODUCTION 

With the growing interest in quantum computation 
and information the search for efficient and robust quan- 
tum gates has become increasingly important. Deutsch 
presented in 1989 a three qubit quantum gate and showed 
that this gate together with arbitrary one-qubit rotations 
are sufficient to create any quantum network [l[. Such 
a set of gates are called universal for quantum compu- 
tation. Since then many sets of gates were proven to be 
universal 0] - one of them consisting of a one-qubit phase 
gate (S), a one-qubit Hadamard gate (H), and a two- 
qubit controlled phase gate (CS). Explicitly in the onc- 
and two-qubit bases ({|0), |1», {|00>, |01), |10), |H)})the 
form of these gates are 
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The purpose of the present work is to show that all these 
gates can be implemented using only gates based on adi- 
abatic evolution using stimulated Raman adiabatic pas- 
sage (STIRAP) and leading to geometric phases. A quan- 
tum system which starts out in the n'th eigenstate of 
a Hamiltonian that changes adiabatically in time will 
according to the adiabatic theorem Q remain in the 
n'th eigenstate but may acquire a phase. In 1984 Berry 
showed that besides a dynamical part 6 n = — J u) n (t)dt 
associated with the eigenfrequency uj n = E n /h the state 
also acquires a geometric phase, 7„ [1] 

*(0) = Vn(0) -> *(t) = exp[i(#„ + 7 n )ty n (t). (2) 

The geometric part depends on the geometric properties 
of the parameter space of the Hamiltonian and Berry 
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showed how to calculate 7„ when the eigenstate is known 
and the system is non-degenerate Q. Berry's phase 
was generalized to degenerate systems by Wilczek and 
Zee [5[ and to non-adiabatic evolution of the system by 
Aharonov and Anandan [||. Quantum computation re- 
lying on these geometric quantum phases is called holo- 
nomic quantum computation 0] and is expected to be 
particularly robust against noise. Previously, proposals 
for holonomic quantum computation were presented for 
nuclear magnetic resonance Q, neutral atoms 0] and 
trapped ions [l(|. 

Control of the geometric phases requires an adiabatic 
evolution and preferably an eigenstate with a zero- valued 
eigenenergy in order to avoid the build-up of an addi- 
tional dynamic phase as the system evolves. A well- 
described adiabatic process is STIRAP where population 
is transferred from one quantum state to another in a 
three-level lambda-system when subject to two different 
laser pulses ordered in a counterintuitive time-sequence 
pd| . The STIRAP process was shown to be very effi- 
cient and robust theoretically as well as experimentally 
[U, [H EE El EE EE EE El- Geometric phases accu- 
mulated during a STIRAP pro cess were previously inves- 
tigated for tripod systems [19( and used for single-qubit 
rotations [2(|, entanglement between atoms in a cavity 
[2l| and holonomic quantum computation with trapped 
ions E^]. In [22| these phases were considered for an 
open quantum system. In this work, we present an ex- 
perimental implementable set of universal gates based on 
geometric phases arising from population transfer in tri- 
pod systems and show explicitly how they depend on the 
experimental parameters. 

The paper is organized as follows: Sec. [Til presents 
the atomic system under consideration. In Sec. IIIII 
we present a one-qubit phase gate and in Sec. IIVI the 
Hadamard gate. The gates are investigated analytically 
as well as numerically. In Sec. [V] we introduce a cou- 
pling between the two qubits and present a two-qubit 
controlled phase gate. In Sec. IVII we discuss the robust- 
ness and conclude. 
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II. TRIPOD SYSTEM 

We consider two atoms (ions or neutrals) with a tri- 
pod level structure as shown in Fig. [T] The three lower 
states (|0),|1) and |2)) are long-lived and coupled to the 
upper state |e) by application of three laser fields with 
Rabi frequencies f^o, ^1, ^2, respectively. In practice 
the lower states can be ground Zeeman- or hypcrfinc- 
sublevels, and |e) is an electronically excited state or ex- 
cited state manifold. We assume the laser fields are on 
two-photon resonance and denote the one-photon detun- 
ing by A. We use STIRAP processes to transfer popula- 
tion among the three lower states and the excited state 
|e) is therefore never populated during the process which 
ensures that no loss occurs due to spontaneous emission. 
We consider {|0), |1)} as our qubit-states. In the follow- 




FIG. 1: Two atomic four-level tripod systems with three laser 
fields applied with Rabi frequencies f2o, Hi, ^2- The one- 
photon detuning is denoted by A, and the subscripts a and 6 
refer to two different atoms. 

ing we implement the universal set of quantum gates ([T]) 
using geometric phases acquired by transferring popula- 
tion adiabatically with STIRAP-processes. 



III. ONE-QUBIT PHASE GATE 

We first consider how to use STIRAP to perform a 
simple one-qubit phase gate: | j) — > e l( ^ j \ j) where \j) is 
either of the two qubit states |0) and |1). For this purpose 
we use a single STIRAP process to transfer the popula- 
tion from \j) to j 2) and another to transfer the popu- 
lation back to \j) again. During this process \j) gains 
a geometric phase. The total pulse sequence is shown 
in Fig. [21 To solve the time evolution during these two 
STIRAP sequences we consider the Hamiltonian for the 
{|j),|e),|2)} lambda system in the rotating wave approx- 
imation (RWA), when we assume two-photon resonance 
between \j) and |2) (see Fig. [T]) 



Ha) = ; 



Ctjit) 
f2*(i) 2 A n 2 (t) 
Q* 2 (t) 



(3) 




FIG. 2: Pulse sequence transferring population from \j) to |2) 
and back to \ j) again. The FWHM of each of the four pulses 
is r (see Eq. (|13[) h AT is the delay between the two sequences 
and At the delay between two pulses within one sequence. 



We parameterize the complex Rabi frequencies as 

% (t) = sin e(t) ^\n 3 {t)\i + \n 2 {t)\\ (4) 
n 2 (t) = coB9(t)^\nM 2 + \Mt)\ i ^ p( - t \ ( 5 ) 

and diagonalize © to obtain the energy eigenvalues 



= A± JA 2 



n 2 



n 



with the non-absorbing zero-valued dark state (uj d 
given by 

\D(t)) = cos 0{t)\j) - sm0(t)e lip(t) \2). 



(6) 



0) 



(7) 



Now we assume that this state is the initial state 
\ip(— 00)) = \D{— 00)) at time t = —00 before the pulses 
and that we vary the real amplitudes of Qj (t) and (t) 
and the phase ip(t) of £l 2 (t) in an adiabatic way such that 
all population stays in \D(t)). During this evolution the 
\D(t)) -state will pick up a phase which is purely geomet- 
ric because the state has zero energy eigenvalue. The 
phase will depend on the evolution of the parameters 9 
and ip and we define these as a vector, R = (0(t), f(t)), in 
parameter-space of the Hamiltonian. Then the acquired 
Berry phase is exactly the integral 

7m =* / ' (D\V R \D) ■ dR = - [ V ^ sin 2 6(t)d<p(t), 

J Ri Jf{ti) 

(8) 

where sin 2 6 is found from Eq. 



sin 2 0(t) 



|Q i (t)| 2 + |Q 2 (t)P 



(9) 



In the pulse sequence of Fig. [5] we assume that all four 
pulses are described by the common function f2(t). In 
addition to f2(i), the Rabi frequency fl 2 (t) is defined by 
the phase <p(t) (see Eq. ([5])). The instants of time t a and 
tb in Fig. [2] are defined such that sin 2 9(t) w for t < t a 
V t > t b + AT and sin 2 6{t) w 1 for t b < t < t a + AT. 
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With these definitions we obtain from Eq. 



In 



ip(t a +AT) 

Idip 

<p(tb) 

V (t b+ AT) n 2 (t + At-AT) 



and © 
(10) 



V (t a +AT) 



fl 2 (t + At- AT) + n 2 (t - AT) 



dp. 



Substituting t' = t— AT in the last integral and assuming 
that ip is a monotonic function we obtain 



7n 



n 2 (t) 



n 2 (t) + n 2 (t + At) 



dip 
~dt 



dt 



(11) 




dt 

n 2 (t' + At) 



n 2 (t' + At) + n 2 {? 
a ^ +AT dip 



■ ^dt' 

dt 



FIG. 3: The upper panel shows the evolution of the popu- 
lation of states \j) (full), \e) (dashed), and |2) (dotted). The 
lower panel shows the evolution of the phases, 4>j of state \j) 
(full) and 4>2 of state |2) (dotted). Numerical and analytical 
results cannot be distinguished on the scale of the figure. The 
calculations were made with sin 2 -pulses (|13|l . <p = t/To and 
The geometric phase thus only depends on the laser field parameters: s2 max ,j /2tt = fi max , 2 /27r = 100/To, At /To = 1, 



dt 



dt = ip(t a )-p{t a +AT). 



phases and requires control of AT, the fraction and 
similarity of the four pulses. All these quantities are rou- 
tinely controlled to high precision in present-day labora- 
tories. After the evolution the final state is 



|j>- 



(12) 



The population and the phases of the three states \j), 
|e) and |2) can be found numerically by solving the time- 
dependent Schrddinger equation. In addition the phases 
can be calculated analytically as described above (Eq. 
(fTTj) and (jT2J)) . In Fig. [3] we show the population of the 
states as well as the evolution of their phases when we 
assume ip = t/To => r y ni = —AT/To and use sin 2 -pulses 
with same amplitude 



Q(t) = 



O max sin 2 ( 




— 1 

2t ' 



if < t < 2r 
otherwise 



(13) 



r/T = 1, AT/To = 5. 



Fig. [2]). The total acquired phase is as shown in Eq. flll[) . 

7m = p(t a ) - p(t a + AT) = -AT/To- The phase of the 
|2)-state (dotted curve) contains not only the geometric 
phase 7 ni but also the additional <p(t) — tt (See Eq. [7]) 
yielding a total phase fa = 7m + <p(t) — 7r. The |2)-state 
therefore accumulates a phase before the first and after 
the second STIRAP process -but non in between where 
7 ni and tp(t) cancel each other. In the time-windows 
where the states in question are populated the direct nu- 
merical solution of the time-dependent Schrddinger equa- 
tion gives results for the phases in agreement with the 
above analytical results. 



The factor of 2 in the pulse assures that r corresponds 
to the FWHM. The populations are shown in the upper 
panel of Fig. [31 Initially all population is in the In- 
state (full curve). During the first STIRAP-process the 
population is transferred from \j) to |2) (dotted curve) 
while the second STIRAP-process transfers the popula- 
tion back to The electronic excited |e) -state (dashed 
curve) is never populated and hence no loss of popula- 
tion occurs due to spontaneous emission from |e). The 
evolution of the phases is shown in the lower panel of 
Fig. [3J The full curve shows the phase of \j). This 
phase is the geometric phase 4>j = 7»i- The phase 
4>j remains zero until the first set of STIRAP pulses 
arrive (t a in Fig. [5]); it then accumulates a phase un- 
til the second pair of pulses has passed (if, + AT in 



IV. HADAMARD GATE 



To implement a Hadamard gate we use three laser 
fields with Rabi frequencies flo, fii and (See Fig.[T]). 
We apply f2o an d f2i in a constant ratio tan#oi = j^j 
with phase difference (/>qi and assume two-photon reso- 
nance. In this case the dressed states of the system are 
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a single dark QDjj)) and two bright states (|+),| — )) 
\D H ) = cos0 O i(t)|O) -sin0 O i(t)e^ ol |l), (14) 
|+) =-L (sm<5(sin0 O i(i)|O) +cos0 O i(i)e l0ol |l)) 
+ cos <5 1 e)) , 

|-) =i (cos<5(sin0 O i(t)|O) + cos0 O i(t)e^»|l)) 
-sin<5|e}) , 



(sin(-|), — cos(-|)) yielding a final state \ipf) = ^(|0) — 

|1)). With this degree of control we therefore produce a 
Hadamard gate with certainty. This gate combined with 
the one-qubit phase gate can generate arbitrary qubit 
rotations and the gate is robust as it depends on con- 
trollable parameters such as the pulse shapes, the ratio 
of Rabi frequencies, the delay between the two sequences 
and the phases of the laser fields. 



with eigenvalues ojd h = 0, uj± = A ± ^ A 2 + + Of 
and 5 defined by tan S = \J • ^ ensure that no 
population is transferred to the |e)-state we choose as 
our initial state 

\i/H) =a\D H ) + V2b(smS\+) + cos5|-)) (15) 



=a\D H ) 
=a\D H ) 



6(sin 0oi | 
b\B), 



COS 0Q1 e 



L ll» 



where a, b are normalization constants and \B) = 
shi6» i|0) + cos6» ie^ 01 |l) is the bright part of \tpi). We 
now use the same pulse sequence as in the case of the one- 
qubit phase gate (Fig. ^ but with both and fii ap- 
plied. The dark state \Dh) does not couple to |2) and is 
therefore unaffected by the pulses, while the bright state 
\B) does couple and is therefore transferred to |2) and 
back again acquiring a geometric phase. In this sense the 
dynamics is similar to the one-qubit case and accordingly 
\B) -> e l "'"H\B) with j riH = - J sin 2 8 H dip H . The geo- 
metric phase 7„ H is controlled by tan0# = ^""q^ 5 -^- 
and ipH, which is the phase difference between fio and 
f^2- After the pulses the system then ends up in 

\4> f ) = a\D H )+be i ^»\B) (16) 
= a\D H ) + 6e i7 "« (sin O1 10) + cos ol e^ 01 11)). 



V. TWO-QUBIT PHASE GATE 



To create a two-qubit phase gate a coupling between 
the two qubits is necessary We consider a coupling 
J5|22)(22|, where E is the coupling strength. In the end 
of this section, we briefly discuss how such a coupling 
can be realized. We assume real Rabi frequencies and all 
laser fields on resonance. If we wish to solve the full sys- 
tem of Fig. [1] analytically it is an advantage to go into the 
interaction picture with respect to Hq = £?|22}(22|. In 
this picture the system has a six-dimensional null-space 
yielding six orthonormal dark states, \Di) (i = 1, . . . , 6): 



\Di) =|00), (18) 
\D 2 ) = -cos0 2 |lO)+sin0 2 |2O), 
|£> 3 > = -cos0 2 |Ol)+sin0 2 |O2), 

\D A ) =-L(sin0 2 (|le) - |el» + cos0 2 (|2e) - |e2»), 

\D 5 ) = cos 2 2 1 11) - sin 2 cos 2 (| 12) + |21» 
+ sin 2 2 e list |22), 

\D 6 ) =^(-sin 2 2 |ll) - sin 2 cos0 2 (|12) + |21)) + \ee) 
v2 

-cos 2 9 2 e lEt \22)). 



To implement the Hadamard gate we control the Rabi 

frequencies to obtain fi max ,2 = y^max.o + ^max.i an< 4 

Oo = — (v2— l)f2i. This latter relation leads to 0oi = J 
and (f>ox = 7r. We control the pulse sequence such that 
7„ H = — 7r and with these parameters we find from (fT4"|) . 
US) and (ig 

\ipi) =(acos|+6sin|)|0) + (asinf -fccosf)|l) (17) 
|V'/) =(acos| -6sin|)|0) + (a sin f + &cos f )|1) 

Now the initial condition \tpi) = |0) corresponds to 
(a, b) = (cos(^), sin(^)) yielding a final state \ipf) = 
4=(|0) + |1)), while \ipi) = |1) corresponds to (a, b) = 



In this degenerate case we use the method described by 
Wilczek and Zee [j| to find the geometric phases. We 
assume that we start with all population in |11) and 
that only n 2 is applied and write ipii - °°) = 
where the index I indicates that we are solving the prob- 
lem in the interaction picture. When we assume an adia- 
batic evolution the population stays within the null-space 
and il>i(t) can be written as (see also Ref. [lj|) 

^(t)=J2B b (t)\D b (t)). (19) 

b 

The time evolution is given by the time-dependent 
Schrodingcr equation (dot denotes differentiation with 
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respect to time) 

]T B b {t)D b {t)) + B b (t)\D b (t)) = 

5 

E B6(t)iA(*)> = - E s f wiaw). (20) 

b 6 

Here we have used that (if(f) - Ho(*))|A>(i)) = for all 
dark states. Taking the inner product with (D c (i)\ yields 

B S) = -J2 B b (t)(D c (t)\D b (t)). (21) 
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The only non-zero (D c (t)|_Dh(t))-elemcnts are 

(D 5 (t)\D 5 (t))=iE sin 4 6, (22) 
(D 5 (t)\D 6 (t)) = -j=E cos 2 9 2 sin 2 2 , 

(D 6 (t)\D 5 (t)) = -^=Ecos 2 6 2 sm 2 6 2 , 

(D 6 (t)\D 6 (t)) = l -Ecos 4 e 2 . 

The differential equations for the ^-coefficients now re- 
duce to 

Bi(t) = 0, B 2 (t) = 0, B 3 (t) = 0, B 4 (t) = 0, (23) 

Bs(t) = -iE sin 4 2 B 5 (t) - -^=E cos 2 6> 2 sin 2 6 2 B 6 {t), 

v2 

i 1 

B 6 (f) = r #cos 2 6> 2 sin 2 2 -B 5 (f) - i-E cos 4 2 B 6 (t). 

v 2 2 

Starting with all initial population in |_D 5 (f)), Eq. |23|) 
shows that |Og(f)) is populated when cos 2 2 sin 2 02 is 
non- vanishing, which is only the case during the turn-on 
and turn-off of pulses. Choosing small pulse widths there- 
fore assures that effectively all population stays in | .05(f)) 
while accumulating the phase. A pulse with FWHM 
t =1 us results in the transfer of 0.5% of the popula- 
tion to |-D 6 (f)) while the smaller FWHM value r =0.5 us 
reduces this population to 0.1%. The value r =0.5 us 
is experimentally feasible but requires higher Rabi fre- 
quencies. Now, in this regime where B$(t) f» 0, we may 
readily solve Eq. |23|) for B§{t) and from [19] we obtain: 

^(f) =e- lf ^-~ sin4e2dt |0 5 (f)) (24) 
= cos 2 e 2 e-' E ^ s ^^ dt \ll) 

-sin 6» 2 cos6» 2 e- 4jE ^ sin4e2dt (|12) + |21» 

■In —iE f* sin 4 2 dt+iEt i o\ 

+ sin u 2 e J -°° 1 22). 



Going back to the Schrodingcr picture an extra phase 
associated with the |22)-energy shift is introduced and 

Ht) =e- l£ l 22 >< 22 lV/(t) (25) 

=6-^/2. -in* fl,* (cos 2 6» 2 |11) 

- sin 2 cos6» 2 (|12) + |21» + sin 2 6> 2 |22>) . 

With all population initially in |11) and application of 
the STIRAP pulse sequence ([5]) to both atoms but with 
only Qi in Fig. [1] applied, we end up in e 7 ™2 111) after the 
pulse sequence, where 7„ 2 = — E J_ f sin 2 dt. If we, 
in addition, keep the phases of the laser fields fixed no 
singlc-qubit phases are accumulated, i.e., 

|00) |00) (26) 
101) -|01) 

|io)-|io) 

|11) -> e i7 "2|n). 

So far, our analysis of the two-qubit gate is general and 
applicable to various atomic systems such as optically 
trapped neutral atoms, trapped ions and rare-earth ions 
doped into crystals. To provide the coupling £'|22)(22| 
we need to be a little more specific, and we suggest in 
the case of trapped atoms or ions to exploit the long- 
range dipole-dipole interaction between Rydberg excited 
atoms [23j and in doped crystals to exploit the interaction 
between excited states with permanent dipole moments 
[2^ | . In [H, [24| the interactions are used for quantum 
gates in stepwise schemes, where first one atom is ex- 
cited, and then the interaction blocks the excitation of 
the second atom, leading to an entanglement between 
them. This stepwise process is not compatible with our 
adiabatic protocol, and we suggest instead to apply the 
interaction to perturb the energy levels of both atoms 
when their |2) states are coupled off- resonantly to the 
excited states. As shown in [25[, this off-resonant exci- 
tation causes an energy shifts (AC Stark shifts) of the 
\2) states, and due to the dipole-dipole interaction, this 
shift will have a non-separable component of precisely 
the desired form for suitable choices of the laser detun- 
ings and strengths. The energy shift E is given by a 
fourth order expansion in Rabi frequencies with a prod- 
uct of three dctunings in the denominator, and the gate 
thus requires relatively long interaction times to avoid 
population transfer to the excited states. 

VI. CONCLUSIONS 

We have shown that in the adiabatic limit popula- 
tion transfer in tripod systems introduces purely geo- 
metric phases. These phases can be used in quantum 
information science to form a set of robust geometric 
gates. The performance of the three gates (flj depends 
on the robustness of the phases: r y ni = — J sin 2 6 dip, 
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ln H = - f** J sin 2 H d(p H and 7„ 2 = -E JJ^ sin 4 9 2 dt. 

The population transfers are done without ever populat- 
ing the upper state |e) in the tripod system (see Fig. [T|), 
which ensures that the gates are insensitive to sponta- 
neous emission. Pulse shapes, delay between sequences 
and ratio between Rabi frequencies are routinely con- 
trolled experimentally without drift in the laboratory. If 
we assume systems where the three laser frequencies lie 
so close that all fields can be generated from the same 
source, the relative phases of the fields are also easily 
controllable. This could be achieved when {|0), |1), |2)} 



are atomic Zccman- or hyperfine-substates. In this case 
all gates are very robust with respect to parameter fluc- 
tuations and can be implemcntablc in present-day labo- 
ratories. 
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